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Abstract 

Recent BABAR measurements on lifetime and mixing of B mesons are reported. Various techniques 
are used, ranging from the full reconstruction of hadronic B decays, to partial reconstruction 
techniques, and to a totally inclusive approach with dilepton events. The results presented are 
based on a data sample collected by BABAR during the 1999-2000 data taking, and should be 
considered as preliminary. 
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1 Introduction 



A precision measurement of the B°B° oscillation frequency is of great importance since this quan- 
tity is sensitive to the CKM matrix element \Vtd\ and, in combination with knowledge of the B®B® 
oscillation frequency, provides a stringent constraint on the Unitarity Triangle. A precision mea- 
surement of B meson lifetimes is also of great importance for the understanding of the dynamics 
involved in heavy quark decays. Moreover, in the specific case of the new asymmetric B-Factories, 
both lifetime and time-dependent mixing measurements are powerful tools for understanding the 
performance of both the detector and the analysis algorithms, and therefore they represent valida- 
tion analyses for the measurements of time-dependent CP asymmetries. 

This paper presents a set of B meson lifetimes and mixing measurements performed with the 
BABAR detector Q at the PEP-II e + e~ asymmetric B-Factory. Several techniques are used, rang- 
ing from fully exclusive to completely inclusive event selections. Exclusive measurements generally 
offer smaller systematic uncertainties with respect to inclusive ones, but they also suffer from smaller 
branching fractions, and consequently larger statistical uncertainties. The following sections will 
detail the three main mixing and lifetime measurement procedures followed in BABAR namely the 
selection of events 

• where one B meson decay is fully reconstructed in a hadronic flavour eigenstate, or 

• containing two high momentum leptons from semileptonic B meson decays (dilepton events), 
or 

• where a B decay is reconstructed by using a semi-exclusive technique in either an hadronic 
or a D*lv final state. 

The data sample used was collected by the BABAR detector at the PEP-II asymmetric B-Factory, 
during the years 1999-2000. Unless otherwise specified, the integrated luminosities are 20.7 fb^ 1 
on the T(4S) peak, and 2.6 40 MeV below resonance. All results are preliminary. 

2 Mixing and Lifetimes with fully reconstructed hadronic B de- 
cays 

This technique is the most similar to the measurement of the time-dependent CP asymmetry also 
presented at this Conference ||, [|, the only difference being the exclusive reconstruction of a 
flavour eigenstate instead of a CP eigenstate. The main ingredients necessary to perform these 
time-dependent measurements are: 

1. to exclusively reconstruct a flavour-eigenstate hadronic final state (B reco ), 

2. to tag the flavour of the other B meson (_Bf a9 )[], 

3. to measure the time difference At between the two B meson decays. 

Since B production at the T(4S) resonance is coherent, tagging the flavour of a B meson at its decay 
time will unambigously determine the flavour of the other B at the same time. The distribution of 

1 This step is not needed for lifetime measurements. 
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the time difference At between B meson decays with the opposite (h + ) or the same (h_) flavour is 
therefore 

h± = l/4r e - r|Ai| (l ± (1 - 2w) cos(Am d At)) (1) 

where T = 1/t b o and w is the probability to get a wrong flavour tagging. The resulting same/opposite 
flavour asymmetry will have a pure cosine time dependence. Equation [j] has to be modified to take 
into account the detector resolution function and backgrounds. 

2.1 Exclusive reconstruction of hadronic B decays 

Several hadronic decays are reconstructed^, based on most favourable production rates and low 
background contamination: 

B° d^+tt- B- -> L>M°7r- 

Dto+p- J/ipK~ 
D(*)+a7/ i){2S)K- 

where D*+ -» D°ir + , D*° -> D°ir°, ip(2S) -» J/ijnc+n- ot£+£~, J/ip -> K*° -► K~tt + . Neu- 

tral D mesons are reconstructed in the K~ir + , K~ir + Tr~ir + , K~7r + ir° and (except for B~ — > D*°tt~) 
KgTT + ir~ modes, charged D mesons in the K~7r + ir + and K^7T + modes. Signal to background dis- 
crimination is based on two variables, the beam energy substituted mass, rriEs, and the energy 
difference AE, defined as 

m ES = \lKlam-VB and AE = E* B - E* beam , 

where E£ eam , p* B and E* B are respectively the beam energy and the momentum and energy of the 
fully reconstructed B meson, computed in the center-of-mass (CMS) frame. The tues resolution 
is about 2.6 MeV and is dominated by the beam energy spread, whereas the resolution on AE is 
mode-dependent and varies between 20 and 30 MeV, being worse in decay chains involving neutral 
pions. The vties distribution for neutral B decays, after applying a 3er cut on AE, is shown in 
Figure |]. The number of neutral (charged) B signal events with tues > 5.27 GeV/c 2 is 6643±96 
(6928±94), with a purity of about 84% (87%). The combinatorial background is estimated by 
fitting the tues distribution between 5.2 and 5.3 GeV/c 2 with the sum of a gaussian for signal and 
an Argus function Q for background. This procedure also allows to assign a signal probability on an 
event-by-event basis, which is used in the fitting procedure. The time dependence of combinatorial 
background is determined from the time distribution of events in the sideband region. In addition 
to combinatorial background, there is a small contribution which peaks in the signal region. This 
peaking background is due to pion/photon swapping between the two B mesons in the event, and 
introduces a contamination from charged B mesons in neutral B decays and vice-versa. Detailed 
Monte-Carlo studies show that peaking background is at the percent level. 

2.2 B-fiavour tagging and particle identification 

After reconstructing a full B meson decay chain (B reco ), the remaining tracks of the event are ana- 
lyzed to determine the Bt ag flavour. Powerful tagging signatures are primary leptons in semileptonic 

2 charge conjugation is always implied in the following. 
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Figure 1: Beam energy substituted mass distribution for fully reconstructed hadronic neutral B 
decays. 



B decays, kaons resulting from b — > c — > s transitions and slow pions from D* in B — > D*X decays. 
The tagging algorithm assigns events to one (and only one) of the following categories, in order of 
decreasing priority: 

• primary lepton tag (electron with p > 1 GeV/c in the CMS or muon with p > 1.1 GeV/c in 
the CMS if no electron is found), with the tag determined by the lepton charge; 

• kaon tag, with tagging defined by the sum of all kaon charges (required to be non-zero); 

• NT1 and NT2 tags, based on the output of a neural network. 

NT1 and NT2 are aimed at recovering leptons which do not pass the requirements for the primary 
lepton tag, resolving ambiguous cases where a lepton tag has a conflicting kaon tag, and identifying 
slow pions from D* decays. Figure |2| shows the tagging efficiency e versus the fraction of wrong 
tags w. The figure of merit for tagging is defined as Q = e(l — 2w) 2 . The kaon tagging category 
has the highest Q, followed by lepton, NT1 and NT2 tags. The number of tagged neutral B events 
is 4538±75, which gives a tagging efficiency of 68%. The mistag fractions are determined from the 
fit to data in the mixing analysis (see 2.4). 



Most of the tagging power relies on particle identification which has therefore to be optimized 
for maximal efficiency and low misidentification probability. 

Electron identification is based on track matching with a calorimeter cluster, on the momentum- 
to-energy ratio (0.89< E/p <1.2, where E and p are the energy measured in the calorimeter and 
the track momentum measured in the drift chamber, respectively), on requirements on the electro- 
magnetic shower shape, and on consistency for dE/dx and the Cherenkov angle measured in the 
DIRC (when available) with the electron hypothesis. 
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Figure 2: Tagging efficiencies vs. wrong tag probability for the four tagging categories described 
in the text. 

Muon identification relies mainly on the Instrumented Flux Return (IFR), namely on the num- 
ber of measured interaction lengths (greater than 2), on the difference between the measured and 
the expected number of interaction lenghts (less than 1), on matching between IFR hits and the 
extrapolated track, and on requirements on the average number and spread of IFR hits per layer. 

Kaons are identified with a neural network based on likelihood ratios computed from quantities 
measured in the DCH and SVT (dE/dx), and in DIRC, where single hits are compared with the 
pattern expected for Cherenkov light in the kaon/pion hypotheses. Kaon/pion separation is above 
3<T for momenta up to 3.5 GeV/c. 

Table |l| shows the efficiencies and pion misidentification rates for electrons, muons and kaons, 
as determined from data control samples. 

2.3 Vertexing and time measurement 

The time difference between the two B meson decays At is inferred from the distance between their 
decay vertices along the beam line, Az, via the boost factor: 

At = Az/ipyc). 

Since < j3"f >=0.56 at PEP-II, the average separation between the two B decay vertices is about 
250 jLtm, which is measurable by a vertex detector. The B reco decay vertex is conceptually easily 
identified, whereas the Bt ag vertex has to be reconstructed inclusively by using all other tracks 
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Table 1: Efficiencies and pion misidentification probabilities for electrons, muons and kaons. 



Particle 


Efficiency (%) 


Pion misid. (%) 


Electrons 


91 


0.13 


Muons 


75 


2.5 


Kaons 


85 
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Table 2: Resolution function parameters for the z position of the fully reconstructed vertex and for 
Az. 



Parameter 


ZBreco 


Az 


Core a(fim) 


45 


100 


Core fraction (%) 


80 


70 


RMS (/xm) 


70 


170 



in the event. Particular care has to be taken in order to avoid any bias coming from tracks not 
originating from the tagging B vertex: 

• neutral long-lived particle (e.g. K°, A) are searched for and used in the vertex fit instead of 
their charged decay products; 

• tracks originating from secondary (charm) decay vertices are excluded from the vertex fit by 
an algorithm which removes iteratively the track which gives the biggest chi-square differ- 
ence between the vertices reconstructed with and without the track itself. The procedure is 
repeated until the chi-square difference for every remaining track is less than 6 or until there 
are no remaining tracks. 

Additional kinematic constraints derived from the B reco momentum and decay vertex position, 
from the beam spot position and size and from the knowledge of the boost, are used during the 
reconstruction of the B tag vertex. 

The resulting Az resolution function is parametrized with three gaussians (core, tail, outlier) for 
the mixing fit, and with a sum of a gaussian and the same gaussian convoluted with an exponential 
decay for the lifetime fit. Table ^ reports the resolution parameters for the z position of the fully 
reconstructed vertex and for Az. The dominant contribution to the Az resolution is due to the 
tagging vertex. In the mixing fit, the widths of the core and tail gaussians are scaled on an event- 
by-event basis with the error computed from the vertex fits, whereas the outlier gaussian has a fixed 
width of 8ps. A possible remaining bias due to secondary vertex tracks not removed in the fitting 
procedure is taken into account by allowing a non-zero mean of the core and tail gaussians. Since 
the bias depends on tagging category, different core biases are allowed for each of them. Most of 
the resolution function parameters are fitted in data: two scale factors for the core and tail widths, 
one tail and four core biases, the relative core/tail/outlier amounts, and a bias for the outliers. 
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2.4 Results 



An unbinned maximum likelihood technique is applied to fit simultaneously the distribution of 
mixed and unmixed events. The signal is parametrized as in Equation whereas the combina- 
torial background is taken into account by a zero-lifetime and a non-oscillatory, non-zero lifetime 
components, both convoluted with a two-gaussian background resolution function. Fit parameters 
are 34: Am^, the mistag fractions uu for the four tagging categories, the signal and background res- 
olution functions, and the background parameters. The correlation between any other parameter 
and Arrid is less than 10%. The Az distributions for mixed and unmixed events and the result- 
ing time-dependent asymmetry are shown in figure 0, together with the fit result. The neutral B 
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Figure 3: Top: At distributions for unmixed and mixed neutral B events. Bottom: time dependent 
asymmetry. The points represent data, the continuous line is the fit result, the dashed line is the 
background contribution. 
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oscillation frequency is measured to be 

Am d = (0.519 ± 0.02CW ± 0.0W syst ) hps' 1 , 

and the systematic uncertainty is dominated by resolution and background modeling. 

Shortly after this Conference, the charged and neutral B meson lifetime measurements presented 
in ICHEP2000 §] were also updated to the full 1999-2000 data sample (figure |). The new results 
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Figure 4: At distributions for neutral (top) and charged (bottom) fully reconstructed B candidates. 
The points refer to data, the line is the fit result, the hatched area represents background. 



arc 



t b0 = (1.546 ±0.032 stat ±0m2 sys t)ps, 
t b+ = (1.673 ±0.032 stat ±0.022 syst )ps, 
t b+ /t b o = 1.082 stat ± 0.026 ±0.011 syst , 



in good agreement with world averages. 
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3 Mixing with dilepton events 



Measuring B meson properties relying only on two leptons from double semileptonic B decays is 
conceptually simpler than a fully exclusive analysis, and more powerful in statistical terms, due 
to the relatively large semileptonic branching fractions and high lepton identification efficiencies. 
Moreover, the B flavour is easily identified by the charge of the lepton. The asymmetry between 
unlike- and like-sign dilepton events evolves as 

N t+ i- (At) - N e+i+/ - e - (At) 

— — : ■ ; = cos AmjAi . 

N en -(At)+N e+e+>e - e -(At) v ; 

However, the main drawbacks of this approach are the presence of non negligible backgrounds due 
to leptons from the b — ► c — ► £ decay chain (cascade decays), which are also the main source of 
wrong tags, and most seriously, the presence of charged and neutral B mesons in almost equal 
proportion in the selected sample, which dilutes the mixing asymmetry. 

Background from cascade decays is minimized with a neural network technique which uses five 
input variables (the lepton momenta and opening angle, and the total and missing energies, all in 
the CMS), and offers better performance in terms of both signal efficiency and mistag rate with 
respect to a traditional approach based on high momentum cuts. 

The fraction R of charged B in the selected sample can be either fit in data or taken (with 
further assumptions) from previous measurements at the T(4S). In the former case, the statistical 
accuracy on the mixing measurement decreases substantially, due to the high correlation between 
the oscillation frequency Am^ and R. In the latter case the systematic uncertainty is dominated by 
the knowledge of R. With the current data sample either approach gives similar total uncertainties. 
However, the present analysis fits R because in the fitting procedure this parameter absorbs any 
difference in efficiency or mistag rate between charged and neutral B decays, which would be 
otherwise another source of systematic uncertainty. 

Other selection criteria of the dilepton analysis include cuts on event shape variables to suppress 
continuum background, on track quality to improve the Az resolution function, and invariant mass 
cuts to reject backgrounds from J/ip decays and photon conversions. The signal purity after all 
cuts is 78%. The main backgrounds are due to events with at least one lepton from cascade decays 
(12%), to events with at least one fake lepton (5%), and to continuum events (5%). 

The z difference between B decay vertices is defined by the z difference of the point of closest 
approach of the leptons to a T(4S) vertex in the transverse plane. This vertex is determined with 
a chi-square fit which uses the two lepton tracks and a beam spot constraint. The resulting Az 
resolution, estimated from signal Monte-Carlo events, has core and tail widths of about 90 fim and 
200 /xm respectively, and 75% of the events in the core. A comparison of the Az distribution for 
leptons from inclusive J/tp decays shows an agreement within 10% between data and simulation. 

The analysis presented in ICHEP2000 || has been updated shortly after this Conference by 
using the entire 1999-2000 data sample. About 100000 events are selected. In the analysis update, 
the total dilepton sample is divided in two independent subsamples, enriched respectively in neutral 
and charged B decays. This is accomplished by reconstructing inclusively the B° — > D* + £~v decay, 
where the D* properties are inferred from the soft pion produced in its decay @. Pions of less 
than 190 MeV/c momentum are identified, the D* energy and momentum are computed and the 
squared missing mass of the neutrino M 2 is computed for both lepton/soft pion combinations. 
An event is assigned to the £>°-enriched sample if a combination satisfies \M 2 \ < 1.5(GeV/c 2 ) 2 , 
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to the i? + -enriched sample otherwise^]. Both subsamples provide similar statistical power to the 
Am d measurement. A simultaneous binned maximum likelihood fit of the same- and opposite-sign 
leptons time distributions on both subsamples gives 

Am d = (0.499 ± 0.010 stat ± 0.012^) hps' 1 , 

which is at present the most precise single measurement of the neutral B meson oscillation frequency. 
Figure |5] shows the time-dependent asymmetries for the two subsamples, together with the fit 
results. The systematic uncertainties are dominated by the knowledge of the resolution function 




Figure 5: At asymmetry between opposite and same sign dilepton events, for the 5 + -enriched (left) 
and i?°-enriched (right) samples. The points represent data, the line is the fit result. 

(0.009), the background parametrization (0.006) and the error on the B meson lifetimes (0.004), 
which are fixed to the PDG2000 || values. Separate fits to the two independent subsamples give 
Am d = (0.498 ±0.015 stai ) hps' 1 and Am d = (0.504 ±0.0U sta t) hps' 1 for the £+-enriched and B°- 
enriched samples, respectively. These results are consistent with an independent analysis, based 
on an unbinned maximum likelihood fitting technique, with the same or even with a different, 
momentum-cut based, event selection. 

4 Overview of other lifetimes and mixing results 

The following results were presented at ICHEP2000 g [n| and are being updated with the full 
1999-2000 BABAR data sample. 

3 if both combinations satisfy the cut, the one with the smallest \M%\ is chosen. 
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4.1 Mixing from fully reconstructed B° — ► _D*£^ decays 

This measurement is very similar to the mixing measurement with fully reconstructed hadronic 
decays, the only difference being in the flavour eigenstate, namely B° — > D*£v with D* — > D°7r + . 
Neutral -D mesons are reconstructed in the X _ 7r + , i , T~7r + 7r~7r + and i^ _ 7r + 7r° decay modes. By 
using a data sample corresponding to a luminosity of 8.9 fb^ 1 at the T(4S) resonance, 7517±104 
B° — ► D*£v are reconstructed and the mixing parameter is measured to be 

Am d = (0.508 ± 0.020 stat ± 0.022 S3/st ) hps' 1 . 

4.2 Lifetimes from partial reconstruction techniques 

Partial reconstruction techniques allow a dramatic enhancement in statistical power, but they are 
more sensitive to the parametrization of physics processes and to biases in the vertex reconstruction, 
the latter due essentially to track swappings between the two B mesons. Two decay chains are 
partially reconstructed: 

• B° — ► £)*7r, where the missing mass due to the non-reconstructed D meson is computed 
by using only the (fast) pion from the B decay, the (slow) pion from the D* decay and the 
appropriate kinematic constraints. About 1700 events are selected by requiring M m i ss >1.854 
GeV/c 2 in 7.4 fb~ l T(4S) data. 

• B° — ► D*£v, where the same slow pion technique already described in Section |3] is applied. 
The signature of a signal event is a pair of oppositely charged lepton-soft pion tracks, whereas 
same sign pairs are used to determine the amount (about 40% of the entire sample) and time 
dependence of the combinatorial background. The analysis selects about 90000 signal events 
on 7.4 fb' 1 T(4S) data. 

The neutral B lifetime is measured to be 

r B o = (1.55 ± 0.05** ±0m syst )ps (B°^D*ir), 
t b o = (1.63 ± 0.02,*,* ±0m syst )ps {B°^D*£u). 

5 Conclusion 

Asymmetric B-Factories offer the unique possibility of precision lifetimes and mixing measurements 
in the B meson sector. Any single measurement performed with the presently available data is com- 
parable to the corresponding world average. The BABAR experiment at PEP-II has investigated 
different reconstruction techniques, ranging from the full reconstruction of a B meson decay chain 
to the totally inclusive selection of events with two leptons from semileptonic decays, to the par- 
tial reconstruction of hadronic and semileptonic decays. Recent preliminary measurements of the 
charged and neutral lifetimes and oscillation frequency with exclusive hadronic decays give 

t b0 = (1.546 ±0.032 stet ±0.022 S2/st )ps, 

r B+ = (1.673 ±0.032 stat ±0M2 syst )ps, 

t b+ /t b0 = 1.082 ±0.026*0* ±0.011^, 

Am d = (0.519 ±0.020 stat ±0.016^) hps' 1 . 
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A preliminary measurement of the mixing parameter with inclusive dilepton events gives 

Am d = (0.499 ± O.OKW ± 0M2 syst ) 

which is at present the most precise single measurement of this quantity. 

Further results on lifetimes and mixing with other techniques are being released in a short time 
scale. 
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